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Abstract 

In this paper we present soliton solutions of two coupled nonlinear Schrodinger equations modulated 
in space and time. The approach allows us to obtain solitons for a large variety of solutions depending 
on the nonlinearity and potential profiles. As examples we show three cases with soliton solutions: a 
solution for the case of a potential changing from repulsive to attractive behavior, and the other two solutions 
corresponding to localized and delocalized nonlinearity terms, respectively. 
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Coupled nonlinear Schrodinger equations (CNLSE) are very important because they can be 
used to model a great variety of physical systems. Solitary waves in these equations are often 
called vector solitons in the literature since they are generically described by two-component wave 
functions [lll,|2|]. One of the simplest vector solitons is known as shape-preserving, self-localized 
solution of coupled nonlinear evolution equations [1]. The CNLSE may also model beam prop- 
agation inside crystals as well as water wave interactions, and in fiber communication system, 
such equations have been shown to govern pulse propagation along orthogonal polarization axes 
in nonlinear optical fibers and in wavelength-division-multiplexed systems [sj]. 

In the case of Bose-Einstein condensates, a many-body description of the effect of Fesh- 
bach resonances in the mean field limit, requires the use of two coupled Gross-Pitaevskii equations 
16,0]. This results in a convenient way to account for the physics of the so-called super-chemistry, 
[81]. In the case of one-dimensional (ID) Bose-Einstein condensates (BECs), the presence of 
cigar-shaped traps introduces tight confinement in two transverse directions, while leaving the 
condensate almost free along the longitudinal axis. The realization of BEC of trapped atoms was 
ly achieved in [9-il|. Next, dark solitons in BEC were formed with repulsive ^^Rb 



experimenta 



atoms in lll2 . 



13n . and bright solitons were generated for attractive ''Li atoms in [Il4l . 



IJ. In the 



case of two-component condensates, the experimental generation has been achieved for different 
hyperfine states in rubidium atoms in a magnetic trap lllql and in sodium atoms in an optical trap 

When treationg a set of CNLSE's, one usually deals with two or more coupled NLSE's with 
constant coefficients. More recently, however, in Refs . lilSMlOll the authors have shown how to 
deal with NLSE with varying cubic coefficient, and in [211] one has extended the procedure to the 
case of cubic and quintic nonlinearities. With this motivation on mind, our goal in the present 
work is to deal with soliton solutions of the CNLSE where the potentials and nonlinearities are 
modulated in space and time. The procedure consists in choosing an Ansatz which transforms 
two CNLSE's into a pair of CNLSE's with constant coefficients, certainly much easier to solve. 
With this procedure we can then obtain analytic solutions for some specific choices of parameters, 
leading us to interesting localized coupled soliton solutions of the bright and/or dark type. We 
illustrate the procedure with several distinct possibilities: firstly, we consider the system with 
attractive potentials and localized nonlinearities. Next, we consider nonlinearities delocalized and 
then we make the potential to change periodically in time from attractive to repulsive behavior. 

The model we start with is described by two CNLSE's, with interactions up to third-order. 
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These equations are given by 



i + + J^(5'ifc(a;,t)|^/'fcnV'i, (1) 



dt dx"^ 

k=l 

I- 



2 



+ t;2(x,t)^2 + 5^('/2fc(x,t)|V'fcnV'2, (2) 



dt dx"^ 

k=l 

where ifjk = tpkix, t), the functions, Vk{x, t), are the trapping potentials, and the Qij, = 1, 2) 
describe the strength of the cubic nonlinearities. Note that we are using standard notation, with 
both the fields and coordinates dimensionless. Note also that we are supposing that our model 
describe ID BEC, and then it should engender cigar-shaped traps to induce tight confinement 
in two transverse directions, leaving the condensate almost free along the x axis. In the super- 
chemistry model of 13,181], the above equations are amended by a term which goes as ^/2aipiip2 
in Eq. ([T]) and with a term of -^^f in Eq. These, number-nonconserving, terms arise from 
treating the Feshbach resonance as a coupling term between atoms and molecules. 

It is interesting to comment here that the elimination of, say, ^2 in favor of an effective equation 
for ipi, would necessarily end up with a GP equation with depletion effect. This could appear in 
the form of complex nonlinearity. We leave such study for a future work. 

In the following we seek to reduce the above two equations into the following pair of equations 

yui^ = + Yl (GikM Ai, (3) 



dC 

d^A 

'dC 



k=l 

d'A, 



2 

/^2A2 = + J2 {G2k\Ak\') A2, (4) 



k=l 

with Hi and the Gij being constants. We achieve this with the use of the twoAnzatse: 

V'i(x,t)=p(x,t)e^^(^'*Mi(C(x,t)), (5) 

^2(x,t)=p(x,t)e^'^(^'*)A2(C(x,t)). (6) 
However, we now have to have 



These equations control the functions p(x, t), ri{x, t) and C{x, t), and the two potentials are given 

by 

We take C{x, t) = F{^{x, t)), with ^(x, t) = x/x{t), which has to obey the Eqs. ©, ® and 
After some algebraic calculations we get that 

9jk = -/^- (11) 

We use Eq. do]) to obtain 



p{x,t) = l/^xF'ia^,t)). (12) 

Also, using Eq. ([8]) we can write 

''^^'^^ = + ^^^^ 

where a = a(t) is a function of time. 

The above results can be used to find finite energy solutions for both ipi and ^2, with appropriate 
choices of F (^) , according to Eq. (fT2)) . We illustrate these results with some applications of current 
interest to BEC. 



Example # 1. - The first example is a extension of the model discussed in feoll . in the case of 
nonlinearities given by 

^?ifc = G,feX"'exp(-3e'), (14) 

with X = x{t) > 0- We consider the simpler case, with pi = fi2 = 0. Here, we get that 
Vi{x,t) = V2{x,t) = v{x,t), with v{x,t) being an attractive potential. We use Eq. (fTTl) to get 
p{x, t) = exp(^^/2) / y/x and the potential is now given by 

v{x,t) = f{t)x^ + h{t), (15) 

with 

1 1 r/^Y 

With the choice a = J x '^dt, we can write h(t) = 0. Thus, the potentials acquire the simple form 
v{x,t) = f{t)x'^. Here we choose Gil = —1/2, G12 = —1,^21 = —1/2, G22 = —1, which we 
use in Eqs. dSHU). In this case, each one of the two solutions, ^jji and ijj2, describes a similar group 
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of atoms with attractive behavior, as is the group of atoms investigated in n2(J\ . With the above 
choice of parameters we obtain the solutions 

' v^dn(AioC, 1/V2) ' 

A.(x,t) = -^-°^"(-^-°^'^/f\ (19) 
2dn(A2oC, 1/v^) 

where Aiq = A20 = Aq = InYJyl j \pl) j ^ , for n = 1,2,... and K{l/y/2) being the elliptic 



integral. The value of Aq is obtained through appropriate boundary conditions; see ll20h . 



We use Eq. (fT2l) to get ( = ^ [1 + erf (^)]. Thus, the pair of coupled solutions are given by 

v^yx dn(A„C,l/^) 

^ ' 2yx dn(^„C, I/V2) 

with being a real function obeying Eq. (fT3]) . 

Next, we use Eq. (fT6l) to obtain the function As two distinct examples, we consider 

the cases / = 1 and / = 1 + ecos((X'ot), to describe periodic and quasiperiodic possibilities, 
respectively. The first case, / = 1, leads us to % = [1 + 15 cos^(2t)]^/^/2. The second case, 
/ = 1 + e cos{uJot), is a little more involved; it leads us to x which is given as a combination of the 
two solutions of the Mathieu equation ^ + 4:f{t)z = 0. It has the form x = (2^1 + 2z|/H^^)^/2, 
where 2 represent the two linearly independent solutions of the Mathieu equation, with W being 
the corresponding Wronskian. 

In Fig.[T]we show the the potential v{x, t) for the two cases, (a) periodic, and (b) quasiperiodic. 
The solution (l20l) is plotted in Fig. |2l for the periodic and quasiperiodic cases, respectively. The 
initial data for the Mathieu equation are z{0) = \/2 and dz/dt{0) = 0. The results for ^2 are 
similar and we do not show them. 

To verify stability of the solutions, we use the split-step method, through finite difference to 
investigate the propagation of the initial states 0) and ^'2(2;, 0) given by Eqs. (I20l) and (|2TI) . 
at t = 0, with the inclusion of a small aleatory perturbation of 3% on the profile of the initial states. 
With this procedure, we could verify that the solutions (|20l) and (|2T]) are both stable against such 
aleatory perturbations. 

Example #2 - In this case, we consider an attractive potential but we allow the nonlinearity to 



change, as we did in Ref. [12 ih . We choose the nonlinearity terms to be, 

9jk = G,kX-'^MiVl^) (22) 



where 7 gives the width of the nonlinearities. With this, we can use Eqs. (fTTI) and (fT2l) to find 

p{x,t) = exp(— ^^/67^)/yx and C(0 = /exp(—,^^/37^)(i^, respectively. The potentials are 
now given by 

V, = f{t)x^ + h{t) - 12, exp(-2eV37'), (23) 

with f{t) e h{t) given by Eqs. (fT6l) and (fTTI) . Here we use a = J x~'^dt, which leads us to h{t) = 0. 
In this case, we can obtain solutions of the form 

= sech(C(x,t)), (24) 

A2(x,t) = -^sech(C(a;,t)), (25) 

with the choice ^1=^2 = Cn = —1, G12 = G21 = —2, and 6*22 = 0, in the Eqs. ((BHU). In 
this case, we note that V'l describes a group of atoms similar to the group considered in [I2 ill . (but 
without quintic interaction) which interacts with the other group of atoms described by ^2- This 
last group, however, does not self-interact. Here, the general solution has the form 

^1(0:, t) = e*^x"'/'exp(-eV67')sech(C(x,t)), (26) 

^2(x,t) = e*''(2x)-'/'exp(-eV67')sech(C(x,t)). (27) 

As before, we see that the temporal evolution of the solutions are directly related to x"^''^' 
which can be obtained from Eq. (fT6l) . For / = 1 and for / = 1 + e cos((X'ot), we can get periodic 
and quasiperiodic solutions, respectively, as show before in the former example. In Fig. [3] we 
show the potential for the case (a) periodic and (b) quasiperiodic, respectively. The solution (l26l) 
is plotted in Fig. HI for the,(a) periodic, and (b) quasiperiodic, cases, respectively. As before, we 
do not depict the case ?/'2, since it is similar to ipi. 

We have studied numerically the propagation of the initial states 0) e ^'2(2;, 0), given by 
Eqs. (I26l) and (1271) . at t = 0, after introducing aleatory perturbation of 3% in the profile of the 
initial states. With this, we could verify that the solutions (|26|) and (l27l) are stable against such 
aleatory perturbations. 

Example # 3 - Next, we consider the case with fii = fi2 = Gn = 0, Gu = —2, G21 = 2, and 
G22 = — 1- We use this in Eqs. dSM]). This leads to two interesting analytical solutions: one field 
leading to a bright-like solution, and the other giving rise to a dark-like solution. They are given 
explicitly by 

= -^tanh(C(x,t)), (28) 



A2(x,t) = sech(C(x,t)). (29) 

We take 

gjk = G,kX-\l + \exp{-e)f, (30) 
with > 0. We use Eq. ^ to get 

p{x,t) = -pr (31) 

(x(l + Aexp(-e)))'/' 

Now, the potentials are given by 

Vj{x,t) = si{t)x'^ + S2{x,t), (32) 

with 

= >^-M-e) A ^ [Aexp(-a-2] A 

This choice of nonlinearity gives Ci^^'t) = ^ + ^/^TXer{{^)/2. With this, the solutions are given by 

^1 (x, t) = -= [x (l + Ae-« j J tanh(C) , (35) 



M^,t) = e*" [x (l + Ae-«')] '^%ech(C), (36) 

with T] real, given by Eq. (fT3l) . and t) as shown above. Due to the several kinds of modu- 
lations of the nonlinearities, we can choose x{i) in several distinct ways. For instance, we may 
choose x(t) = 1 + a sin(t) + /3 sm{\/2t), with a and (3 being constants related to the periodic or 
quasiperiodic choice of x- 

In Fig. ©, we plot the potential (l32l) for the (a) periodic and (b) quasiperiodic cases. Since these 
potentials have some structure at small distances, in Fig. Q we show the same potentials, but now 
for small x. The parameters were chosen to be a = 0.1 and /3 = for the periodic solution, and 
a = /3 = 0.1 for the quasiperiodic solution, with A = 0.5. Despite the presence of structures 
at small distances, the main importance of the potentials are their long distance behavior, which 
allows the presence of localized excitations, even though the potential changes from attractive to 
repulsive behavior periodically. We get to this result after changing A — t- — A, which significantly 
alters the small distance behavior of the potential, but preserving its long distance behavior; even 
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in this case, the system supports similar soliton solutions, showing that the small distance behav- 
ior of the potential does not affect the presence of localized solutions. This is clearly seen in 
Figs. dV]) and dS]) where we show the corresponding solutions and |?/'2p, for the periodic and 
quasiperiodic choices of considered above. 

Due to the presence of a dark-like solution, which does not vanish in the limit x — t- ±oo, the 
numerical method used in the former two examples fails to work in such a case. This, therefore, 
requires seeking alternative numerical algorithm to circumvent this problem. We are currently 
investigating this issue. 

In summary, in this paper we have presented soliton solutions for two coupled nonlinear 
Schrodinger equations modulated in the space and time. This model is robust in the sense of 
presenting a vast quantity of nontrivial solutions in the systems with modulations in space and 
time of the traps and nonlinearities. We have considered three different examples of potentials and 
nonlinearities. All these solutions are of great interest in several areas, such that, in coupled Bose- 
Einstein condensates and communications in optical fibers. Our work can potentially motivate 
future studies and help guide possible experimental work in superchemistry and fluid dynamics in 
general. 
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FIGURE CAPTIONS: 

FIG. 1 : Plot of the potential vi for the cases for (a) periodic and (b) quasiperiodic, for e = 0.5 

e coq = 1. 

FIG. 2: Plot of iV^ip (Eq. (I20l) ) for the two cases (a) periodic and (b) quasiperiodic, for n = 1, 

e = 0.5 and luq = 1. 

FIG. 3: Plot of the potential f (x, t) (Eq. (|23l) ) for the cases (a) periodic and (b) quasiperiodic, 
for 7 = 6, e = 0.5 and uq = 1. 

FIG. 4: Plot of l^/'ip (Eq. (1261) ) for the (a) periodic and (b) quasiperiodic, for 7 = 6, e = 0.5 
andcjo = 1- 

FIG. 5: Plot of the potential v{x, t) (Eq. (fTOl) ) for the cases (a) periodic and (b) quasiperiodic. 

FIG. 6: Plot of the potentials show in Fig. [5] for x small. In (a) we display the case periodic, 
and in (b) the quasiperiodic case is shown. 

FIG. 7: Plot of l^/'ip (Eq. (1351) ) for the cases (a) periodic and (b) quasiperiodic. 

FIG. 8: Plot of |^/'2p (Eq. (|36l) ) for the cases (a) periodic and (b) quasiperiodic. 
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FIG. 1: Plot of the potential vi for the cases for (a) periodic and (b) quasiperiodic, for e = 0.5 e = 1. 
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FIG. 2: Plot of iV'iP (Eq. (l20l )) for the two cases (a) periodic and (b) quasiperiodic, for n = 1, e = 0.5 and 

UJQ = 1. 
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FIG. 3: Plot of the potential v{x,t) (Eq. (|23]) ) for the cases (a) periodic and (b) quasiperiodic, for 7 = 6, 
e = 0.5 and coq = 1. 



13 




0.4 0.8 1.2 1.6 2 0.5 1 1.5 2 2.5 



FIG. 4: Plot of IV'iP (Eq. (l26l) ) for the (a) periodic and (b) quasiperiodic, for 7 = 6, e = 0.5 and ujq = 1. 
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FIG. 5: Plot of the potential v{x, t) (Eq. (ITOl )) for the cases (a) periodic and (b) quasiperiodic. 
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FIG. 6: Plot of the potentials show in Fig. |5] for x small. In (a) we display the case periodic, and in (b) the 
quasiperiodic case is shown. 
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